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Twelve years ago, Haldane formulated his famous conjecture for 1-d antiferromagnetic quantum spin chains. 
In the context of the 2-d 0(3) a model with a term, it predicts a phase transition as —> 7r, which has not yet 
been verified reliably. To simulate this we use the Wolff cluster algorithm together with an improved estimator for 
the charge distribution. Each cluster carries integer or half integer charge. Clusters with charge 1/2 are identified 
with merons. At = ir they are inactive, such that the mass gap vanishes. We obtain critical exponents which are 
consistent with predictions from the k = 1 WZNW model, therefore confirming a second order phase transition. 



1. Haldane's Conjecture and Numerical 
Results 

Antiferromagnetic quantum spin chains are 
most easily described by the Heisenberg Hamil- 
tonian 



H 



(1) 



where J > and (ij) are nearest neighbors. We 
deal with quantum spins, therefore 
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This implies that the Neel state, \s, — s, s...), is 
not an energy eigenstate (for classical spin chains 
it is the ground state). In fact, the low energy 
behavior is rather complicated. 

In 1983, Haldane conjectured that there is a 
qualitative difference between integer and half in- 
teger s [1]. For half integer spin he predicts a 
gapless phase, which would have been expected in 
general from the spin wave theory. The surpris- 
ing point was that for integer s, he conjectured a 
finite mass gap. 

For spin chains, Haldane's conjecture has been 
proved analytically for all half integer spins [2] 
and for s = 1 [3]. The latter case also agrees 
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with neutron scattering experiments, which found 
a mass gap in the quasi 1-d s = 1 antiferromag- 
net CsNiCi3 [4], in contrast to s = 1/2 systems 
such as CuCl2- 2 (NC5H5). Moreover there is nu- 
merical evidence for s = 2 [5], but not yet for 
2<s£ where the assumed mass gap is tiny. 

A low energy effective description (or a large s 
expansion) leads from the 1-d quantum spin chain 
to the 2-d 0(3) non-linear a model with a 6 vac- 
uum angle [1,6] 



6 = 2tt(s mod 1). 



(3) 



Here the poorly understood case flips: while the 
mass gap at 6 = is well known, its vanishing 
at 6 = 7r has no strong evidence so far. A nu- 
merical study by Bhanot et al. did not arrive at 
conclusive results [7]. Haldane and Affleck have 
proposed a meron picture to describe this situa- 
tion. They advocated a second order phase tran- 
sition at 6 = 7r. Affleck extracted quantitative 
predictions for the critical exponents by an iden- 
tification with the k = 1 WZNW model. 

The action of the 2-d 0(3) non-linear a model 
in Euclidean space reads 

S[e\ = ((fx^d^e-d^e-^-e^e^d^ex d v e)](A) 

where |e*(a;)| = 1 and by large s construction from 
the antiferromagnetic quantum spin chain one ob- 
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tains 6 = 2-ks. The term (I/871-) J (Pxe^e- (d^ex 
d v e) counts how many times e(x) covers S 2 . It is 
the topological charge Q £ Z . Writing the action 
as 

S[e\ = S [e\-iOQ[e\ (So £ IR) (5) 

we see that physics is 2w periodic in 6, so we may 
restrict 6 e.g. to [0, 2w[ and obtain (3). 

Haldane and Affleck have introduced the fol- 
lowing picture of the mass gap at 6 = due 
to meron condensation [8]. Include a potential 
(m 2 /2g) (e 2 ) 2 . Thus e tends to the xy plane. A 
vortex in this plane has Q = 1/2 and can there- 
fore be called a meron. If the distance r from 
the origin is infinite, then the vortex should lie 
strictly in the xy plane, in order to achieve fi- 
nite action. Close to the origin, however, it has 
to leave this plane to avoid UV infinite action, 
in particular e z (r = 0) = ±1. Thus merons and 
antimerons do provide finite contributions to the 
action and can disorder the system, i.e. cause a 
mass gap. At m = 00 we have the 0(2) model 
with the Kosterlitz-Thouless phase, understood 
by vortex condensation. As we decrease m down 
to 0, the massive phase should cover the whole 
range of the coupling g. 

Now we include the topological term. A 
meron/antimeron picks up the factor exp(±i#/2). 
In particular at 6 = w they cancel; they become 
inactive and can no longer cause a mass gap. 

This picture, however, is far from rigorous. In- 
troducing m and sending it to at the end is a 
dangerous interchange of limits. If we let m = 
from the beginning then the vortex picture does 
not make any sense at all. Moreover there are 
multi-vortices which do not cancel at 6 = w. 

In the region addressed by this picture, (3 = 1/g 
larger than some minimum, we expect a second 
order phase transition. Since we refer to solid 
state physics, there is no need for scaling in the 
continuum limit (3 — > 00. Conformal field theory 
renders the above picture more quantitative [9]. 

The k = 1 WZNW theory describes a class of 
models with a massless phase, in particular an- 
tiferromagnetic half integer quantum spin chains, 
the Hubbard model and the 0(3) model at 6 = w. 
The central charge k must be an integer for the 
path integral to be well-defined, k = 1 is the only 



case free of relevant operators, therefore it is likely 
to be a stable fixed point for the above models and 
it predicts their critical exponents. Marginal op- 
erators add logarithmic corrections, which further 
improve the agreement with our data, see below 
and [10]. In our case we obtain a mass scale 

m{6) oc |tt - 6\ 2 ' 3 ■ \ ln{\w - 6>|) I" 1 / 2 . (6) 

Now we can express the topological and magnetic 
susceptibility Xt, Xm in terms of Fisher's dimen- 
sionless parameter z = Lm (L: size of the system) 

Xt{6) = (L/V^L) g t (z), (7) 
Xm {0) = WlnT g m {z). (8) 

Close to the phase transition, g t and g m are sup- 
posed to be universal functions. Our data clearly 
confirm this. Note that a first order phase tran- 
sition is ruled out in the regime we simulate, 
since in that case the susceptibilities would di- 
verge oc L . Fig. 1 shows g%, g m for volumes 
V = 18 x 6, 24 x 8, 30 x 10 and 36 x 12. We 
simulated at extremely strong coupling, (3 = 0, 
such that the correlation length at 6 = is only 
2.8. For each volume we performed at least 10 7 
sweeps. In the plots the measured points are in- 
terpolated, without any fitting of free parameters. 

2. The Algorithm and its Interpretation 

The simulation of models with a 6 vacuum term 
is notoriously difficult. Since a complex Boltz- 
mann factor is not a probability, ordinary Monte 
Carlo fails. In addition one usually has a sign 
problem as 6 approaches w. 

We simulate at 6 = and then incorporate 
the factor exp(i#<3). This requires an extremely 
precise determination of the charge distribution 
p(Q), which yields the partition function Z{6) 
[11], 

p(Q)=[De8 Q , Q[el e- s °W, Z{6) = Y J V{Q)e ieQ '(9) 
J Q 

Observables at fixed charge and 6 = 0, and finally 
expectation values at any 6, are given by 

(V)Q = J De6 Q , Q[el O[e\e- s °W, 
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Figure 1. Universal functions fift(-z) and g m (z) 
based on numerical data. 
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For technical reasons we deal with a triangu- 
lar lattice with periodic boundary conditions im- 
posed by a hexagon of side length L, such that 
there are 3L 2 (independent) sites. We attach to 
each site a unit vector e x . The charge of one tri- 
angle is the minimal solid angle spanned by the 
vectors at its corners. The total charge Q £ Z is 
the sum over all triangle charges. 

We also specify the lattice action at 6 = 0: 

S [e\ = ^ s{e x ,e y ) 
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c 2/ > 2 

otherwise 
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i.e. the permitted relative angles of nearest neigh- 
bors are < 2w/3. We are in the universality class 
of the it model, as we recognize from the sym- 
metries. A priori this lattice action is as good as 



the standard action. The reason for imposing the 
constraint is purely technical, see below. 

We apply the Wolff cluster algorithm [12]. We 
first choose a random unit vector r. 'Flipping' a 
spin e x — > e x ' then means to reflect it in the plane 
through perpendicular to r. Now we connect 
neighboring spins e x , e y by a bond with proba- 
bility 





s(e x ,e y ) — s(e 



») 



£y ^ &x 

otherwise 



These bonds form clusters. The spins in one clus- 
ter can only be flipped collectively, but the clus- 
ters are independent. If we flip one cluster, the 
total charge changes Q — > Q' . We define the 
charge of this cluster as Q c = (Q — Q')/2, an in- 
teger or half integer. Now our constraint has the 
highly non trivial consequence that the charge of 
one cluster is independent of the orientation of 
all other clusters. In a typical diagnostics (for 
V = 36 x 12, g = oo) most clusters are neutral 
and have a mean size of (n c ) = 1.5 sites; about 4 
percent have \Q C \ = 1/2, (n c ) = 10; 0.1 percent 
have \Q C \ = 1, (i c ) = 32 and very few clusters 
have higher charges. 

The fact that the cluster charge is determined 
locally enables the construction of an improved 
estimator for p(Q) and M 2 (Q) (the magnetiza- 
tion squared), which are needed for Xt an d Xm, 
respectively. A configuration with N clusters is 
related to 2^ configurations by cluster flips. Sim- 
ple combinatorics yields the charge distribution 
for the whole ensemble. Typically we have O(100) 
clusters, hence the gain in statistics is tremen- 
dous. Furthermore Z(w) only gets positive con- 
tributions from each ensemble, so there is no sign 
problem in this algorithm. For Xm we have to fix 
Q, but still we can flip all the neutral clusters. 
Thanks to the improved estimator, the improve- 
ment factor for the statistics is still huge. 

p(Q) varies enormously, but we need to know it 
precisely over many orders of magnitude. Here we 
use a trial distribution pt(Q) (as close as possible 
to p(Q)) together with a reweighting technique 
[11]. We work with an effective action 



S efl [e\ = S [e\ + In p t {Q[e\) 



(12) 
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related to the charge distribution by 

p(Q)=Pt(Q) J De8 Q , Q[g] e- s '«W. (13) 

A good pt(Q) makes the charge distribution in 
S e f[ approximately flat. In general it is hard to 
find, since one does not know p(Q) in advance. 
But in our case the improved estimator permits 
the determination of p(Q) to a good accuracy over 
25 orders of magnitude, as shown in fig. 2. 
Pt{Q) is then used for flipping the clusters ac- 
cording to the Metropolis algorithm. 
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Figure 2. Topological charge distribution 

This cluster formulation not only provides a 
very efficient algorithm, but also suggests a topo- 
logical interpretation. Let us consider the clus- 
ters as physical objects, the crucial degrees of 
freedom. Then clusters with Q c = 1/2 are in- 
terpreted as merons. Since the clusters are in- 
dependent, at 6 = the merons build an ideal 
gas, which causes the mass gap. Now we include 
the factor e i9Q ; at 6 = ir it is (-1) Q . If a con- 
figuration has at least one cluster of half integer 
charge, then flipping that cluster changes the sign 
factor, so the ensemble containing this configura- 
tion does not contribute to Z(w). The merons 
and antimerons are 'bound in pairs', they can no 
longer disorder the system and the mass gap van- 
ishes. 

This is a precise formulation of the picture of 
Haldane and Affleck. The advantage is that it 
works without breaking the 0(3) symmetry. 



3. Outlook 

Seiberg has shown that at very strong coupling 
there is a phase transition of first order at 6 = w 
[13]. We simulated at 1/g = /3 = with the 
constraint imposed in (11), which corresponds to 
a P s > in the standard lattice action. Since 
we still observed second order, we could run our 
algorithm at (3 < and try to reach the first order 
regime. A new study claims agreement with first 
order, but is plagued by large errors at 6 w [14]. 

It is straightforward to apply our algorithm to 
the Ising model. This project is in progress. 

Finally the topological interpretation of the 
clusters might even give us a hint for a suitable 
definition of clusters in lattice gauge theory. 

We thank I. Affleck for helpful comments. 
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